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Executive  Summary 


The  subject  of  tne  research  reported  here  is  the  numerical  simulation  of  hypersonic 
turbulent  flow  from  first  principles  utilizing  the  full  time-dependent  compressible  flow 
Navier-Stokes  equations.  The  objective  of  this  Phase  I  SBIR  project  was  to  demonstrate 
the  feasibility  of  using  such  a  mathematical  formulation  to  expose  the  essential  physics 
of  turbulence  unique  to  hypersonic  flow.  Motivation  for  proposing  this  project  was  derived 
from  intense  interest  by  the  U.S.  Air  Force  in  the  TAV  and  NASP  programs,  as  well  as 
on-going  work  in  hypervelocity  projectiles.  Research  which  is  currently  being  carried  out 
by  The  Hokenson  Company  in  this  last  area  for  AFATL/MNG,  Eglin  AFB  has  guided  the 
configuration  of  the  fiowfield  which  was  simulated  during  Phase  1. 

The  viscous  and  heat  conducting  flow  of  a  perfect  gas  along  right  circular  cylinder  whose 
axis  was  aligned  with  the  undisturbed  upstream  flow  was  computed  over  a  range  of  free- 
stream  Mach  numbers  from  i.4  through  6.3.  The  cylinder  could  be  rotated  about  its  axis 
so  as  to  establisn  a  mean  velocity  field  with  three  components,  albeit  axisymmetric. 
Without  rotation,  the  boundary  layer  mean  flow  was  two-dimensional  in  the  conventional 
sense.  Tne  upstream  boundary  condition  was  specified  in  a  unique  way  so  as  not  to  intro¬ 
duce  excessive  numerical  'noise'  into  the  simulation.  The  technique  developed  herein 
was  to  allow  the  (axial  and  azimuthal)  flow  to  'transition'  from  full  velocity  slip  at  the 
wail  to  zero  slip  over  a  finite  number  of  axial  grid  points.  This  allowed  the  boundary 
layer  to  develop  naturally  without  generating  the  huge  disturbances  that  are  computed 
if  the  upstream  uniform  flow  abruptly  encounters  a  zero  slip  surface.  This  tactic  also 
significantly  reduced  the  time  and  costs  associated  with  solving  the  entire  upstream  flow 
(bow  shock,  stagnation  region,  transonic  region,  etc.)  associated  with  a  particular  fore¬ 
body  geometry,  which  was  not  the  focus  of  this  brief  research  effort. 


In  this  manner,  supersonic  and  hypersonic  fiowfieids  which  were  exceptionally  quiet 
were  computed.  For  even  very  large  Reynolds  numbers,  the  initial  numerical  noise  level 
was  so  low  that  the  spatial-temporal  instability  growth  envelope  forced  unacceptably 
long  computational  times.  Therefore,  the  computed  boundary  layer  flow  was  excited 
artificially  with  an  harmonic  disturbance  on  the  surface  of  the  cylinder  at  several  grid 
points  well  upstream  of  the  observation  location.  Both  (temporally)  periodic  heating  and 
cooling  as  well  as  suction  and  blowing  were  studied,  the  latter  being  adopted  for  presen¬ 
tation  here  over  the  range  of  Mach  numbers  of  interest.  During  Phase  II  various  ambient 
noise  models  shall  be  added  to  the  simulation  for  more  natural  excitation. 

Numerical  simulations  of  the  entire  flowfield  were  carried  out  for  a  single  value  of  the 
wall  mass  transfer  disturbance  magnitude  (equal  to  1%  of  the  freestream  mass  flux)  and 
various  frequencies  relative  to  that  associated  with  the  viscous  flow  length  and  time 
scales.  The  results  indicate  a  distinct  Mach  number  effect  when  the  azimuthal  vorticity, 
entropy  and  pressure  fluctuations  are  observed.  In  the  supersonic  flow  regime  the  viscous 
and  inviscid  disturbance  fiowfieids  are  fully-coupled,  with  the  entropy  fluctuations  pro¬ 
pagating  well  beyond  the  location  where  the  vorticity  fluctuations  terminate.  This  is 
reminiscent  of  the  wavelets  observed  experimentally  emanating  from  TBL  structures. 

At  higher  Mach  numbers,  the  entropy  fluctuations  also  extend  beyond  the  vorticity  fluc¬ 
tuations  but  do  not  encompass  the  entire  inviscid  flow.  Instead,  a  multiple  layer  picture 
appears  in  which  two  boundary  layer  edges  may  be  defined  relative  to  the  vorticity  and 
entropy  fluctuations. 

As  a  result  of  this  work,  the  feasibility  of  using  first  principle  numerical  simulations 
of  hypersonic  turbulence  to  expose  the  unique  physics  of  a  boundary  layer  flow  has  been 
demonstrated.  In  addition,  it  is  likely  that  the  results  of  such  computations  couid  also 


be  useful  for  developing  more  compact  engineering  representations  of  the  turbulent 
transport  processes.  Based  on  intense  Air  Force  interest  in  the  aforementional  projectile 
and  aerospace  vehicle  projects,  it  is  clear  that  an  extensive  computational  study  carried 
out  in  a  6.1  Phase  II  SBIR  effort  would  produce  information  and  data  that  could  be  tran¬ 
sitioned  rapidly  for  use  in  various  6.2  programs  and  beyond. 


Discussion 


Introduction 


During  the  past  few  years  the  resurgence  of  research  activity  in  the  field  of  hypersonics 
has  been  motivated  by  U.S.  Air  Force  interest  in  the  TAV  and  NASP  program,  as  weJJ 
>  as  various  high  speed  projectile/missiJe  applications.  Based  on  extensive  experimental 

studies,  the  space  shuttle  was  developed  more  than  a  decade  ago  without  a  satisfactory 
comprehensive  theoretical  formulation  for  the  turbulent  transport  processes.  Instead,  a 
patchwork  auilt  of  semi-empirical  analyses  were  tuned  to  the  available  data  in  order  to 
predict  full-scale  heat  transfer,  drag,  separation,  etc.  Unfortunately,  in  the  intervening 
years  between  the  shuttle  R&D  and  the  current  hypersonic  flow  programs,  breakthroughs 
\  in  the  understanding  of  turbulent  shear  flow  structure  have  so  far  been  focused  in  the 

lower  speed  flow  regimes.  As  a  consequence,  the  accurate  prediction  of  heat  transfer 
rates,  for  example,  in  turbulent  hypersonic  flows  is  still  essentially  an  art  form  guided 
)  by  the  available  data.  Numerical  predictions  of  the  time-averaged  flow  is  carried  out 

within  the  context  of  modeling  principles  known  to  be  generally  valid  only  in  quite  diffe¬ 
rent  Mach  number  regimes.  Relative  to  the  knowledge  of  turbulent  structure  physics  in 
,  lower  speed  flows,  where  transport  models  can  be  relatively  successful,  our  understanding 

of  the  detailed  processes  in  hypersonic  shear  flows  is  indeed  minimal. 

,  For  incompressible  shear  flows,  numerical  simulations  are  now  being  used  to  predict  the 

detailed  vortical  dynamics  in  turbulent  structures.  This  is  important  not  only  in  and  of 
itself  to  understand  nature  and  have  the  ability  to  predict  chaotic  flows  from  first  prin- 
.  ciples,  but  also  in  providing  a  basis  for  developing  representations  of  various  fluctuation 

correlations  in  time-averaged  flow  transport  models.  In  supersonic  flow,  experimental 


data  regarding  the  detailed  temporally-  and  spatially-resolved  structure  of  turbulent 
vorticity,  entropy  and  pressure  fluctuations  is  orders  of  magnitude  less  than  that  avail¬ 
able  in  incompressible  flow.  In  addition,  first  principle  computation  of  such  phenomena 
is  in  its  infancy,  whereas  the  prediction  of  time-averaged  flow  using  various  semi-empiri¬ 
cal  models  has  been  in  high  gear  for  decades.  The  success  which  such  an  approach  has 
enjoyed  rests  on  the  abundance  of  time-averaged  data  with  which  models  can  be  'tuned' 
as  well  as  the  apparently  simple  structure  of  at  least  the  small  amplitude  disturbances, 
as  first  exposed  by  Kovasznay  and  his  co-workers: 
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The  utility  of  dilitation  as  a  variable,  in  lieu  of  entropy,  has  been  exploited  by  the 
Principal  Investigator  in  previous  research.  The  enclosed  paper  indicates  that  the  struc¬ 
ture  of  non-linear  fluctuations,  such  as  those  encountered  in  hypersonic  flow,  may  be 
more  readily  understood  when  dilatation  is  selected. 
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Dilatation  wave/soliton  dynamics  in  variable  density  flows 

Gustave  Hokenson 

Dynunucs  Tecnnftcgy.  Inc..  Torrent <?.  (V  Jorma  9C50I 
i Received  24  July  1981;  accepted  for  publication  5  January  1982* 

The  dynamics  of  dilatation  has  been  investigated  for  variable  density  flows  •  -hich  exhibit  a 
temporal  evolutionary  character  In  lieu  of  decomposing  the  conservation  of  mass  equation 
(appropriate  tor  strictly  nor.diffusive/ineompressiblc  flows;  and  utilizing  div  u  =  0  to  generate  a 
Poissur  pressure  equation,  the  closure  relation  “small  div  u"  was  used.  This  suggests  an 
approximate  decomposition  of  the  div  u  governing  equation  into  a  steady  elliptic  pressure 
equation  and  a  time-dependent  equation  for  nonacoustic  waves  which  governs  dr-  u  therein 
denoted  by/}  ..  to/?  /or.  —  [1  -  \Sld.x  -  1 1  /  4  ReilO  lr:  p/dx>\.dfl  fdx\  —  O  '  —  i  l/ReiicF/2  /dxr\. 
utilizing  the  exact  equations  of  continuity  and  momentum.  The  solution  of  this  j"v  orticity-like"i 
equation  prov  ides  information  on  how  an  initially  small  div  u  field  evolves,  frequently  into  a 
solitary  wav  e.  Should  the  dilatation  remain  “small."  rather  than  be  inonphysically  forced  equal 
to  zero,  the  full  conservation  relations  may  be  utilized  However,  if  div  u  grows  large,  the  elliptic- 
pressure  equation  is  inconsistent  with  the  untampered  governing  equations.  The  motivation  for 
this  approach  is  to  provide  a  structuie  vv  ithm  which: !  1 1  small  diffusive/compressibility  effects 
may  be  studied.  <2-  the  significance  of  div  u  physics,  in  flows  which  are  traditionally  treated  as 
approximately  dilatation  free,  may  be  analy  zed,  and  :?s  evolutionary  pressure  transport 
formalisms  may  bt  developed  and  evaluated  for  computational  efficiency/accuracy  and  the 
inclusion  of  nonlinear  acoustic  waves. 

PACS  numbers  4'  iO.  -  g.  4“.3G  —  s.  47  55. Hd 


INTRODUCTION 

Fluid  flows  which  are  approximately  mcempiessible 
and  nondifiusive  are  often  analyzed  computed  utilizing  a 
decomposed  conservation  of  mass  equation  Dp  DtzzO  and, 
therefore.  /2~ 0  [see  Lq  Mi.  Appendix]  This  allows  the 
governing  system  of  equations  to  be  closed,  independent  of 
dealing  with  the  equations  of  stale  and  thermodynamics. 

In  this  context,  incompressible  means  that  the  equation 
of  stale,  e  g  .  p  -  p\f.  5  ,  is  independent  of  pressure  and 
relaie-  the  fluid  density. p.  to  the  local  thermodynamic  tem¬ 
perature.  T.  and  various  chemical  species  constituents.  S 
1  he  term  nondiffusivc  indicates  that,  relative  to  macroscop¬ 
ic  flowfieid  length  scales,  the  scales  characteristic  of  molecu¬ 
lar  transport,  over  typical  flow  times,  are  negligible  and, 
therefore.  DT S  Wr  Clearly,  that  implies 
Dp  Dt re- 0  and,  from  Eq  iA  1 1.  0. 

The  arguments  required,  to  establishes,  'Dt^O  involve 
only  small  transport  coefficients  and  zero  species  production 
via  chemical  reaction  However.  DT  /DrcsO  only  if  heat  and 
momentum  transport  coefficients  are  small  and  it  is  assumed 
a pnon  that  buoyancy  and  pressure  work  is  negligible 
A  consistent  Dp  D: z-  0 f>  foundation  is  thus  established 

A  complementary  way  of  viewing  this  traditional  equa¬ 
tion  decorr.pc.s'.tior,,  which  relates  to  subsequent  discussions, 
involves  differentiating  Eq.  ,AI..  resulting  m 

DS ! 

—  =  n  -  -  k.  :  1 1 

Di 

w  her ipK  ~  D  y:Dr  According  Eq  I if  k  =  0  exactly, 
by  virtue  of  the  af- .'emeriti.  :.ed  assumptions.  SI  -  {>  initial 
condition-  would  be  maintained  for  all  time  However,  jn 
arb;;rariU  small  initial  dilatation  anomaly  >11,  approacnes 
v.  as  Uj—\.  in  conve.ted  coordinates  A  strictly  in- 


eompressible/nondifiTusivc  model  would  require 
H  =  /?,.  =  0  as  a  requisite  auxiliary  relationship  on  the  basis 
of  boundedness  alone.  However,  imposing  this  constraint  on 
slightly  diffusive/compressible  flows  has  a  significant  effect 
on  the  mathematical  structure  and  the  physics  which  it  re¬ 
flects.  as  indicated  by  Eq.  i  1  i. 

In  the  context  of  Dp;  Dt^  0,  utilization  of  the  full  con¬ 
servation  relations  to  solve  for  an  evolutionary  J?  field  is. 
therefore,  apparently  a  delicate  matter  which  is  the  focus  of 
this  work.  The  sensitivity  of  its  mathematical  structure  and 
Eq  'll  itself  are  reminiscent  of  the  vorticity  transport/ 
stretching  equation  In  fact,  the  symbol  Q  was  chosen  to 
represent  do  u  precisely  to  emphasize  the  i  ~  1  character  as  o> 
conventionally  does  for  vorticity  Icurl  u!.  For  homogeneous 
flows,  vorticity  is  typically  generated  locally,  e  g.,  at  solid 
boundaries,  shock  wave  intersections,  etc  ,  and  diffuses/'con- 
vects  into  adjacent,  nearly  irrotational  regions  Similarly, 
the  div  u  field  equation  serves  the  variable  density  flow  by 
representing  the  generation  of  dilatation  in  certain  diffusive ' 
compressible  "boundary  lay  ers"  in  space/time  and  its  subse¬ 
quent  influence  on  neighboring,  nearly  dilatation-free  flow. 
Such  "boundary  layers"  might  arise  due  to  the  presence  of 
large  gradients  near  sources  of  heat  and  chemical  species. 
The  associated  vorticity  production  term  in  nonhomogen- 
eous  flows.  Xp  >  Vp  '.  is  also  active  in  these  dilatation 
boundary  layers 

EQUATION  DEVELOPMENT 

In  order  to  expose  the  "vorticity-like"  role  which  div  u 
may  play  in  both  the  phy  sics  and  computational  treatment  of 
rotational,  variable  density  flows,  the  following  formulation 
was  developed  It  is  intended  initially,  to  illuminate  the 
unique  features  of  dilaiational  dynamics  and,  subsequently . 
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to  provide  a  framework  within  which  the  accuracy  and  gen¬ 
erality  of  the  analy  sis  may  be  increased 

Combining  the  governing  equations  of  motion1  (Al) 
and  i A2|.  subsequent  differentiation  and  summing  leads  to 
the  following  exact  expression  for/2: 


a!  incompressible/nondiffusive  model  t/2  =  Oi  of  variable 
density  flows  reduces  Eq  (2 1  to 


d'-p  td\np\  dp 

dx,dx  v  ax  )  dx, 

d2u,uJ  _  fd  In  p  \  d2u, 

^  dx  dx  \  dx,  /  dx.dx, 


13) 


Note  that  setting  fl  —  0,  in  lieu  of  a  small  value,  has  dramatic 
consequences  on  the  mathematical  structure  of  Eq.  I2i,  con¬ 
vening  an  evolutionary  equation  into  an  elliptic  equation  in 
space  This  structure  is  consistent  with  the  physical  interpre¬ 
tation  of  infinite  signal  propagation  speed  in  dilatation-free 
flows.  In  addition.  Eq  1 3 1  displays  explicit  flow  nonuniform- 
ity  effects  ld\np/dx,\  which  modify  the  classical  Poisson 
pressure  equation  by  introducing  pressure  gradient  and  vis¬ 
cous  terms,  “weighted"  by  the  density  variation. 

In  lieu  of  dealing  with  a  precise,  vet  inexact,  closure 
relation  !/?  —  Oi  and  decomposing  conservation  of  mass  into 
approximate  equations,  it  is  of  interest  to  investigate  an  alter¬ 
nate.  imprecise,  closure  technique  in  the  context  of  the  exact 
conservation  of  mass  relation.  As  will  be  shown,  the  effect  of  a 
"small  /2  "  closure  is  tv.  suggest  the  decomposition  of  Eq.  (2) 
instead  of  i A 1  i.  The  stability  and  faithfulness  to  the  "small 
/2  "  assumption  of  this  approach  is  the  subject  of  this  study. 

Realizing  that  fl  =  0  is  often  a  relatively  weak  assump¬ 
tion  in  variable  density  flows,  the  sum  of  all  P  terms  in  Eq.  \2) 
is  estimated  to  be  of  order  e,  where  f<l  u/L  |;  and  \u/L  i,  1  is 
a  characteristic  macroscopic  fiowfield  time  scale.  The  feasi¬ 
bility  of  retaining  Eq. |3|.  along  with  the  untampered  govern¬ 
ing  conservation  equation  of  mass  and  momentum  is 
investigated 

Note  that  the  essence  of  this  argument  is  identical  to 
that  previously  used  to  decompose  the  conservation  of  mass 
equation  according  to  D  In  p/Dt  —  c  and  P  =  —  f.  where 
e  =  0  is  implemented  after  evaluating  the  smallness  of 
D  In  p  *  Dt  from  the  equations  of  state,  thermodynamics,  and 
species  In  the  present  case,  each  group  of  terms,  namely 
those  /?  terms  m  Eq  i 2 .  and  the  pressure/velocity  terms  in 
Eq  :3i,  is  smali  and  subsequently  set  equal  to  zero  The  con¬ 
sequence  of  this  closure  approximation  is  to  establish  gov¬ 
erning  equations  [Eqs  'A  1 1,  i  A2  .  '3i]  w  hich  implicit >  force 
/2  to  vary  according  to 


dt  V  3  dx,  }  dx,  3  dx , dx 


To  make  this  development  more  precise,  various  formal  ex¬ 
pansions  of  Eq.  i2i  are  possible  However,  for  the  purposes  of 
this  work,  the  following  heuristic  approach  is  sufficient 
Consider  Eq  (2)  in  tne  representational  form 
<9/2  /dt  -rj\f)\  =  e,  w  here  e  is  the  error  arising  from  the  fact 
that  Eq.  I3|  is  not  satisfied  exactly  [in  Eq  I2i],  It  is  assumed 
that  f  is  small  and  /2  in  Eq.  ;2i  is  expanded  in  the  form 
fl  =  22,  -r  f/2;  -t-  H  O  T.  To  first  order.  Eq  i4;  is  retained 
which  essentially  filters  out  all  nonlinear  acoustic  aspects  of 
dilatation  wave  propagation 

To  reiterate,  Eq.  (4i  is  a  consequence  of  solving  the  un¬ 
tampered  conservation  of  mass  and  momentum  equations  in 
conjunction  with  an  elliptic  pressure  field  relation  In  gener¬ 
al,  Eq.  :4|  need  not  be  solved  since  its  only  function  is  to  serve 
as  an  fl  growth  monitor.  However,  the  /2  characteristics 
may  be  assessed,  to  a  considerable  extent  independent  of  spe¬ 
cific  fiowfield  details,  by  analyzing  the  solutions  of  Eq.  |4|. 

Due  to  nonzero  initial/boundary  conditions,  fl  will  ex¬ 
hibit  some  evolutionary  behavior.  If  it  maintains  a  “relative¬ 
ly  small"  value,  the  solution  of  the  governing  equations  will 
remain  faithful  to  the  untampered  conservation  relations. 
However,  if  fl  grows  inonphysically),  the  formulation  for 
pressure  expressed  in  Eq.  |3|  is  inconsistent  w  ith  the  small  fl 
closure  model,  at  least  in  the  context  of  the  full  conservation 
relations,  and  an  alternate  pressure  formulation  is  required. 


ANALYSIS 

In  order  to  examine  the  characteristic  behavior  of  Eq. 
(4),  consider  a  one  spatial  dimension  case  for  which  u—>ua  as 
x~*  _  x  _  The  ensuing  expression  may  be  nondimensiona- 
lized  according  to  Eqs.  |A3),  resulting  in 


2£_(,+ 

dr  V  )  4  Re  Sx  1  dx 


.  1 

-n  !S| 

Initially,  it  is  of  interest  to  consider  the  solutions  to  Eq  (5)  for 
various  limiting  cases.  Note  that  a  uniform  22  field  vanes 
according  to 


and,  as  with  Eq.  II),  approaches  as  as  22„r— *1.  Similarly,  a 
nonuniform  22  field  evolves  according  to  Eq.  (6)  as  Re—*  oc  in 
characteristic  coordinates  defined  by 

dx/dr  =  -  f  1  +  f Qdx - —  |7) 

V  J  4  Re  dx  J 


At  least  for  these  limits,  the  lineanzed  solutions  moti¬ 
vated  by  situations  w here  220<  1  degenerate  into  solutions  of 
the  form  fl  =  220.  Considering  the  fact  that  arbitrarily  small 
values  of  22(l  become  unbounded  for  sufficiently  large  r,  these 
solutions  are  of  tnvial  importance  in  the  assessment  of  the 
limitations  of  the  "small  22  "  closure  model  proposed  here 
For  cases  with  finite  Reynolds  number,  the  linearized 
version  of  Eq.  ( 5 1  becomes 

(922  (92?  1  d:Jl  ,0 

dr  dx  Re  dx’ 
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where  —  (1  —  d  In  p/d.i/A  Re)  In  this  case  the  fi  solu¬ 
tion  takes  the  form  /}  „  c'  Kc '  ‘  ‘  ;  so  that  the  solution 

simply  diffuses  relative  to  the  curve .v  —  cr/2  =  0.  However, 
sufficiently  large  Reynolds  numbers  invalidate  the  viscous 
linearization,  even  for  small  with  the  behavior  indicated 
by  Eq.  (61  taking  ov  er  in  either  the  intrinsic  or  characteristic 
coordinate  system.  Therefore,  simultaneous  consideration 
of  both  viscous  and  nonlinear  effects  is  s’enerally  required 
As  an  illustrative  example  which  has  some  bearing  on 
subsequent  computations,  consider  the  case  for  which  the  fi  ' 
nonlinearity  is  dominated  by  \\Sldx\\dfi  /d.ii  In  this  situa¬ 
tion.  Eq.  1 8:  becomes 

dfi  -  dfi  1  d'fi 

- rCifit  —  = - — r- .  |9| 

dr  d.x  Re  dx~ 

where  c=  —  ll  —  ffidx  —  d  lnp/<?x/4  Rei 

It  is  of  interest  to  search  for  steads  solutions'  in  the 
convected  frame  of  reference  defined  by  ?/=-.x  +  i'r.  where 
i '  is  a  iconstanti  wave  convection  speed  In  this  coordinate 
system  Eq.  (9;  becomes 

ici/?  ■  f  \n  =  —  /?-.  dOi 

Re 

where  the  prime  denotes  differentiation  with  respect  to  r/. 
Equation  !  10.  may  be  integrated  to  provide 


the  primesi  approach  zero  as  i  and  i)—  +  ac . 

Equation  1 14)  may  be  integrated  twice  to  prov  ide 

J—^-L~if=A.T}  +  A„  115) 

Re  2 

where  /.4  and  /.«  are  constants  of  integration.  Applying  the 
boundary  conditions  on /to  Eq.  Il5l  results  in 

—  +  — I  [/=  1+  t,  (16) 

Re  2  * 

which  may  be  integrated  in  the  form 

f - e - — - -  Re  (dr),  (1?) 

J  -f-/2±Lf+ll/2^L't  J 

resulting  in 

/=  ±  C— £  tanhl  —  <5  Re  77/2),  (18) 


/?  =  <5  :lRe/2)sech:f  —  <5  Re  77/2),  (19) 

w  here 

6=IU:^2U -r  l),  :.  (20) 

Note  that,  since  sech  is  symmetric  about  77  =  0  with  a  maxi¬ 
mum  value  of  1.0,  Eq.  (19)  provides  the  relationship 

=6  2  Re/2,  (211 


Cl  (Hi 

-  _  Re 

where  C  =  jcdfi  and  is  an  integration  constant.  Various 
"jump-type"  phenomena  may  be  exposed  by  seeking  expo¬ 
nential  solutions  (approximately  1  consistent  with  the  func¬ 
tional  dependence:  Sfidx  =  -  (2/.2)tf.  where  A:  is  a  con¬ 
stant  In  this  case,  C  may  be  written.  (/  ,fi  -t-  /.  J?  :  •+-  /.  ,1, 
w  here  A  ,==  —  1  -r  d  \np/dx/\  Re  and  Eq.  ( 1 1|  may  be  inte¬ 
grated  in  the  form 

Re  77  =  f  — - 7 - — — = - 7 — ,  (12) 

J  (Xn  —  /.  ,)  -e  (A ,  ±  U 1/2  +  A2fi ' 

to  provide  the  solution 

n.  —  n  ~ 

-J- — =  exp  (/.;  Re  fflflj  —  fi,))-  (13) 

It  is  clear  from  Eq.  (13)  that,  as  77—*  ±  K,fi—fi,f  exponen¬ 
tially.  Therefore,  the  “jump”  in  fi  across  this  solitary  wave  is 
fir  -  fi,  and  fi . ,^are  related  to  the  parameters  of  the  prob¬ 
lem  by  fi.fi,  =  [S.n  —  A-,)/a:  and  fi,  1-  fif 

-  —  U|  ±  V i/A,  Additionally,  as  77—*  ±  x.Eq  (10) may 
beevaluated to provide  +  U  =  [C 1  fif)  —  C \fi,  )]/lfij  —  fi  ). 
Therefore,  for  certain  situations,  the  presence  of  "shock"  or 
“jump-type”  solitary  wave  solutions  to  Eq.  (5i  should  be 
anticipated 

In  order  to  assess  the  fully  nonlinear/viscous  case,  con¬ 
sider  Eq  1 5:  with  the  generally  small)  density  nonuniformity 
term  temporarily  set  equal  to  zero.  Once  again,  hypothesiz¬ 
ing  similarity  solutions  of  the  form  fi  =  fi\-q\.  where 
-7  Crand  T'isa  iconstanti  wave  convection  speed.  Eq 
5i  may  be  written 

r  if  -JT  -  -^-r  <'4' 

Re 

The  dimensionless  velocity,  /ixi,  is  defined  such  that 
df  idx  =  /}  and /-—l  while  all  its  derivatives  las  indicated  by 


and,  therefore 

t/=2i±l±l2/?ma>/Re|,':|,  (22) 

where  the  first  ±  sign  corresponds  to  the  one  in  Eq  (14), 
whereas  the  second  arises  from  the  two  roots  to  the  solution 
At  this  point,  either  root  may  be  physically  reasonable  with 
the  "observed'  value  dependent  on  the  details  of  the  initial 
conditions  and  its  evolution.  However,  it  is  clear  that  nega¬ 
tive-valued  initial  conditions  are  incompatible  with  this  soli¬ 
tary  wave  solution.  The  <5  =  0  solutions  admit  negative  val¬ 
ues  of  fi  but  require  fi— * «  as  77— 0.  Note  that  constant 
density  nonuniformity  cases  may  be  accommodated  by  rede¬ 
fining  U  to  be  U  +  d  In  p/d i/4  Re. 

Therefore,  for  certain  ranges  of  Re,  arbitrarily  small  fin 
values  grow  nonlinear  and  form  various  types  of  stationary 
(in  the  convected  reference  frame)  waves,  which  remain  a 
permanent  part  of  the  solution. 

COMPUTATIONS 

As  discussed  in  ’he  previous  section,  given  any/20,  there 
exists  a  Re  such  that,  for  sufficiently  large  r,  the  fi  solution 
becomes  dominantly  nonlinear.  In  lieu  of  utilizing  small  val¬ 
ues  of  fi„  and  computing  long  evolutionary  periods,  the  ap¬ 
proach  taken  here  is  to  expose  the  unique  features  of  dilata¬ 
tion  dynamics  by  focusing  on  larger  fi„  values  with 
consequently  smaller  critical  Re  and  nonlinear  evolution 
times  Therefore,  Eq.  (5)  has  been  integrated'  for  various 
Rey  nolds  numbers  and  initial  conditions,  assuming  a  Gaus¬ 
sian  fi,,  disturbance  at  a  given  spatial  location  The  function¬ 
al  form  used,  exp  —  4i:|,  reflects  a  characteristic  length 
which  equals  the  width  of  the  initial  conditions  at  the  e- 
folding  amplitude.  Analysis  indicated  that  the  nonuniform¬ 
ity  term  Id  In  p/dx/ 4  Re)  is  of  secondary  importance  (rela¬ 
tive  to  the  dominant  terms,  1  +  Sfidx]  and  has,  therefore, 
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FIG.  1  The  evolution  of  a  Gaussian  initial  O  disturbance  of  magnitude  /?„ 
and  (d|  Re  »  4.0. 


:  1.0  at  a  Reynolds  number.  |a)  Re  «  1.0.  (b|  Re  »  2.0.  (c|  Re  < 


been  set  equal  to  zero  in  these  calculations.  The  results  of  the 
integrations  illustrate  that  both  bounded  and  unbounded  O 
solutions  arise. 

The  effects  of  Reynolds  number,  initial  condition  am¬ 
plitude  and  sign  are  depicted  in  Figs.  1-4.  For  a  positive 
initial  O  disturbance  of  amplitude  1 .0,  the  effect  of  Reynolds 
number  is  depicted  in  Figs.  1(a)— 1  (d).  Note  that  the  dilatation 
disturbance  wave  propagates  upwind.  Therefore,  O  will  dis¬ 
tort  the  field  ahead  of  its  origin  and  eventually  contaminate 
the  entire  domain.  At  Reynolds  numbers  of  approximately 
10.0  or  smaller,  the  initial  amplitude  diminishes  with  time 
until  an  approximately  steady  profile  is  attained  in  the  con- 
vected  frame  of  reference  at  r~3.0,  as  discussed  in  the  pre¬ 
vious  section.  For  Reynolds  numbers  greater  than  10.0,  the 
initial  conditions  grow  with  time. 

At  a  Reynolds  number  of  one,  the  effect  of  an  increas¬ 
ingly  positive  initial  condition  was  computed  and  is  present¬ 


ed  in  Figs.  2(a)-2(d).  Note  that,  with  increasing  amplitude, 
both  nonlinearities  in  Eq.  (5)  come  into  play  with  increasing 
strength.  This  is  evident  from  the  increased  amplitude  (rela¬ 
tive  to  O0)  and  its  upwind  “tilt”  The  solution  exhibits  a 
decrease  in  amplitude  with  time  for  initial  amplitudes  less 
than  ~  10.0.  As  the  amplitude  increases,  the  convection 
speed  increases  and  tilts  the  solution  to  the  left,  as  the  ten¬ 
dency  for  a  “growth-mode”  solution  is  approached  at  an 
amplitude  of  approximately  10.0 

When  the  sign  of  the  initial  conditions  is  changed,  die 
effect  of  the  nonlinear! ties  on  the  solution  is  altered-  As 
shown  in  Figs.  3(a)— 3(d)  for  an  amplitude  of  —  1.0,  increas¬ 
ing  the  Reynolds  number  narrows  the  spatial  extant  of  the 
disturbance  at  large  times  with  a  dramatic  change  in  its  char¬ 
acteristic  shape.  Note  that  the  initially  rapid  damping  of  the 
initial  conditions,  doe  to  the  (now  negative)  fi1  terms,  de¬ 
creases  and  the  disturbance  eventually  relaxes  toward  aaro 
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FIG  2  Evolution  of  >  Gaimian  initial  /3  disturbance  of  lal  mafnitude  f5„.  -  2.0  at  a  Reynolds  number.  Re  -  1.0:  |b)  magnitude  /?«. .  -  3  0  at  a 
Reynolds  number.  Re  -  10.  (cl  magnitude  -  4  0  at  a  Reynolds  number.  Re  -  1.0;  (d|  magnitude  ■  5.0  at  a  Reynolds  number.  Re  «  10. 


asymptotically.  In  this  case,  the  effect  of  a  change  in  the  sign 
tof  fftdx  is  to  direct  the  wave  downstream  and  to  eventually 
dominate  the  nonlmearities  by  modifying  the  shape  accord¬ 
ing  to  whether  <9/7  /dx  is  a  source  or  sink  term  in  the  overall 
dynamic  balance. 

Finally,  for  a  Reynolds  number  of  10.0  (above  which 
the  solution  changes  only  slightly)  the  effect  of  increasingly 
negative  initial  conditions  is  presented  in  Figs.  4(s)-4(d).  As 
the  initial  conditions  become  more  negative  several  conse¬ 
quences  on  the  solution  are  observed.  First,  although  the 
initial  decay  of  the  solution  is  increasingly  effective  in  damp¬ 
ing  the  disturbance,  its  spatial  extent  is  enlarged.  In  addition, 
the  characteristic  shape  becomes  a  more  nearly  flat  profile 
with  relatively  abrupt  transition  to  the  undisturbed  frees- 
stream  at  its  leading  and  trailing  edge,  remisisoeat  of  the 
“jump-type”  solitary  wave  profile  discussed  in  the  previous 
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section.  Finally,  the  abruptness  at  the  leading  edge  of  the 
wave  propagating  downstream  is  significantly  different  than 
the  relatively  smooth  trailing  edge  waveform. 

CONCLUSION 

It  has  been  shown  that,  when  the  untampered  conserva¬ 
tion  equations  of  mass  and  momentum  are  solved  in  con¬ 
junction  with  a  steady  elliptic  pressure  field  equation,  com¬ 
plicated  dilatation  dynamics  arise  which  significantly  alter 
the  computed  flowfield  characteristics.  Depending  on  the 
Reynolds  number,  sign  and  amplitude  of  the  initial  condi¬ 
tions.  both  bounded  and  unbounded  solutions  exist  for  the 
dilatation  evolution.  For  positive  initial  conditions,  an  up¬ 
stream-propagating  solitary  wave  solution  exists  which  may 
contaminate  the  entire  flowfield.  However,  its  amplitude  is  a 
fraction  of  the  disturbance  initial  condition.  For  negative 
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FIG.  3.  Evolution  of  a  Gaussian  initial  O  disturbance  of  magnitude  fJ„ 
and  (dl  Re  =  10.0 


=  -  1.0  at  a  Reynolds  number,  (a)  Re  *  1.0,  (b|  Re  -  2.0,  |c)  Re  «  3.0, 


initial  conditions,  the  initial  disturbance  spreads  out  and 
forms  two  “jump-type"  transitions  connected  by  a  “roof¬ 
top"  profile.  Therefore,  the  Poisson  pressure  equation  is 
compatible  with  the  full  conservation  relations  only  within  a 
window  determined  by  the  parameters  which  characterize 
the  dilatation  disturbance.  Outside  of  this  window,  a  "small 
/?  "  formulation  utilizing  the  untampered  conservation  equa¬ 
tions  is  inconsistent  with  Eq.  (3),  due  primarily  to  the  invali¬ 
dation  of  the  arguments  which  led  to  the  decomposition  of 
Eq.  (2|. 

On-going  work  involves  alternate  forms  for  an  evolu¬ 
tionary  pressure  (transport)  equation  which  complements 
Eq.  (2j.  The  objective  is  to  alleviate  both  the  steady  pressure 
equation  inconsistency  and  the  need  for  generating  a  compu¬ 
tationally  time  consuming  (yet  approximate)  elliptic  field  so¬ 
lution,  thus  admitting  nonlinear  acoustic  waves.  Along  with 
the  generation  of  an  elliptic  field  pressure  equation  goes  the 


establishment  of  a  set  of  equations  which  exceeds  the  num¬ 
ber  unknowns,  albeit  not  linearly  independent.  Therefore, 
some  machinations  are  required  to  satisfy  all  equations  si¬ 
multaneously,  typically  via  an  ft  “correction"  term  to  Eq. 
(3).  Inasmuch  as  the  “small  (l "  decomposition  of  Eq.  (2) 
avoids  the  equation  excess  problem,  it  is  felt  that  increasing 
the  accuracy  of  an  /2 -oriented  approach  may  be  fruitful  in 
generating  a  completely  evolutionary  pressure  formalism  for 
nonacoustic  as  well  as  acoustic  waves.  Short  of  embarking 
on  an  exact  variable  p  computation,  the  conservation  of  mass 
equation  suggests  the  following  form  for  a  pressure  transport 
relation: 

+  pfl  /e'  —  0,  (23) 

where  f'mw-1  for  isentropic  flow  and  a  is  the  local  speed  of 
sound.  For  gases,  in  this  case  p/2  /('  m  ypQ,  whereas  in  iso- 
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FIG  4  Evolution  of  t  Gtussian  initial  O  disturbance  of  |a|  magnitude  /2_.  -  -  2.0  at  a  Reynolds  number.  Re  *  10.0:  (b|  magnitude 
»  -  3.0  at  a  Reynolds  number.  Re  =  10.0,  (cl  magnitude  »=  —  4.0  at  a  Reynolds  number,  Re  —  10.0;  (d)  magnitude  “  —  5.0  at  a  Reyn¬ 

olds  number.  Re  =  10.0. 


thermal  (or  DT  /Z>f~0)  flow  pfl  /('  =  p/2.  For  nondiffusive 
(T and  S, )  liquid  flows  e"  =  {dp/dp)sj  which,  in  general,  is  a 
function  of  the  local  S,  and  T,  forcing  theS,  and  T  fields  to  be 
computed.  In  tome  cases,  however,  dp /dp  for  a  liquid  may  be 
approximately  a  function  of  p  alone. 

Regardless  of  the  particular  f'  formulation,  the  cou¬ 
pling  of  Eqs.  (23)  and  (21  forms  a  candidate  evolutionary 
pressure  transport  system  to  mimic  the  predictions  of  Eq.  ( 3), 
while  admitting  nonlinear  acoustic  waves,  and  allows  the 
conservation  of  mass  equation  to  be  retained  in  a  computa¬ 
tionally  effective  framework  The  unique  (l  dynamics  indi¬ 
cated  in  Eq.  (2)  would  be  similar  to  those  encountered  here  in 
Eq.  (5|  and,  therefore,  the  present  computations  serve  as  a 
model  for  the  behavior  of  Eqs.  (2)  and  (23).  An  appropriate 
form  for  e'  must  be  validated.  However,  for  moderate  /2,  this 
structure  might  allow  Eq.  (3)  to  be  replaoed  without  t1  re¬ 
flecting  the  detailed  compressibility/dilfusive  physics  In  the 


framework  of  exact,  multiequatioo  solutions  for/2,  the  effect 
of  coupling  through  -  (1/4  Re)d  lap/dx  is  readily  examined. 
However,  the  structure  of  the  solitary  waves  is  not  then  ana¬ 
lytically  accessible,  due  to  the  form  of  -  (1/4  Re)  (ln/5)': 
(ft  /4  ReKl  +  S  fidx  ±  V),  derived  from  continuity.  At  Re 
of  interest  this  effect  is  generally  small. 

In  addition,  the  implications  of  complex  12  dynamics 
[not  only  in  Eqs.  (4),  (3),  and  (6),  but  also  in  Eq.  (23)  where 
small  i 2  is  amplified  by  e]  to  problems  of  physical  impor¬ 
tance,  which  are  traditionally  treated  within  the  (2  «  0 
framework,  is  being  studied,  e.g.,  in  vortex  stretching  dyna¬ 
mics.  Finally,  utilization  of  this  technique  to  approximately 
quantify  the  significance  of  coeapressibihty /diffusive  phe¬ 
nomena  in  weakly  nonhomogeneous  Sows  is  being  investi¬ 
gated.  With  the  inclusion  of  noaBnaar  acoustic  waves,  the 
effect  of  bulk  viscosity  on  the  dflatation  phenomena  in  liq¬ 
uids  and  nonmonatomic  gases  must  also  be  considered.  In 
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Note  that,  for  linearized  supersonic  flow,  there  are  no  first-order  interactions  between 
the  vorticity,  entropy  and  pressure  fluctuations.  In  hypersonic  shear  flows  the  available 
turbulent  structural  data  is  essentially  nil  and  it  is  likely  that  the  intrinsic  nature  of 
the  vorticity,  entropy  and  pressure  fluctuations  is  more  complex  than  in  supersonic  flow. 
Unfortunately,  the  interest  in  and  need  to  predict  hypersonic  turbulent  flows  is  far  out¬ 
pacing  our  ability  to  acquire  data  or  develop  adequate  compact  mathematical  formula¬ 
tions  thereof. 

Therefore,  in  spite  of  the  fact  that  numerical  simulation  of  turbulent  structure  in  even 
low  supersonic  Mach  number  shear  flows  is  still  an  evolving  art,  an  attempt  to  predict 
hypersonic  turbulent  vorticity,  entropy  and  pressure  fluctuation  dynamics  is  indicated 
on  several  bases.  First,  a  general  first  principle  prediction  investigation  will  provide 
fundamental  new  information  on  the  essential  flowfield  structures/physics  which  are 
unique  to  this  regime.  At  the  same  time,  such  simulations  shall  serve  to  guide  forthcom¬ 
ing  experiments  on  hypersonic  turbulent  transport  processes.  Equally  important,  however, 
is  the  role  which  such  simulations  could  play  in  the  formulation  of  compact  mathematical 
representations  of  various  fluctuation  correlations  required  in  time-averaged  flow  trans¬ 
port  mooels  of  TAv,  NASP  and  various  hypervelocity  projectile/missile  flowfields.  The 
pressure  of  U.S.  Air  Force  interest  in  such  programs  will  not  allow  hypersonic  turbulent 
flow  R&D  to  await  the  natural  evolutionary  acquisition  of  experimental  data,  develop¬ 
ment  of  first-principie  prediction  of  hypersonic  vorticity/entropy/pressure  fluctuations, 
and  the  eventual  construction  of  compact  transport  models.  Instead,  the  aforementioned 
R&D  programs,  an  a  others,  require  the  timely  execution  of  a  high  payoff  research  pro- 
gra^  which  attempts  to  leao-frog  some  of  the  traditional  steps  of  evolution.  The  Phase 
I  SB1P  effort  reported  upon  here  was  proposed  in  order  to  demonstrate  the  feasibility 
of  using  direct  numerical  simulations  of  hypersonic  shear  flows  as  a  vehicle  to  assist  in 
this  leap-frogging  ard  thus  achieve  the  aforementioned  goals. 


Formulation 


The  numerical  simulations  carried  out  during  this  Phase  I  SBIR  effort  employed  the 
full  time-dependent,  viscous  and  compressible  flow  Navier-Stokes  equations  written 
in  the  following  vectorized  form: 


-17- 


|1.  I'a  1^1  |l.  .  i '  .  (>.  1>,  ll.  t  ,  il.  |l.  |<.  i'aj'l  i1.  !<a  (  J 


A4  ■.«  »■»  I 


The  Hokenson  Company 


EQUATIONS  OF  MOTION 


dU  bF  bC  bH 
—  +  —  +  •—  +  —  -0 
bi  bx  by  dz 


L'=  •  pi 


puu-  T,, 


fivi-TU 


i  P  i 


,  pEu  -  A  —  -  ( «:  , ,  4-  it  , , 
ox 


*  wrU  ) 


/3 M’,  bWt  2  aM,  \ 
rv-h'>P~,l\Jx--Jxi  Vv  bX>  ) 


j  PUW-T,, 


PI'*  —  T;t 


PWH'-Tjj 


pEl  -X—  -(  UT,;  +  VT!}  +  HT,;  )  ; 


p£n  -A  —  -  (  U7  j ji  +  I’T.V  +  H'Tji 


GUSTAVE  J  HOKENSON.  PhD*  CHIEF  SCIENTIST*  840  S.  TREMAINE  AVE.,  LOS  ANGELES.  CA.  90005  *TEL.  (213)  935-3743 

-18- 


Oil 


S 


K 


B 


I 


E 


fi 

feS 


I 


s 

S' 


r 

% 

$ 


p; 


a 


•s: 


The  gas  was  assumed  to  follow  the  perfect  gas  law  and  obey  the  Newtonian  and  Fourier 
formulations  for  molecular  transport  of  momentum  and  heat.  Due  to  time  constraints, 
no  extensive  competition  amongst  the  various  available  numerical  schemes  in  the  exis¬ 
ting  code  was  carried  out  in  order  to  determine  the  optimum  approach  for  this  problem. 
As  a  result,  MacCormack's  explicit  scheme  was  selected  and,  on  the  basis  of  various 
test  solutions,  the  traditional  artificial  viscosity  terms  were  deleted.  An  exponential 
compression  of  the  grid  spacing  near  the  wall  allowed  the  shear  flow  to  be  fully  resolved. 

In  the  course  of  evaluating  the  aforementioned  numerics,  various  geometries  were  inves¬ 
tigated.  Due  to  the  problems  which  were  encountered  with  planar  geometries  (of  both 
finite  and  infinite  span-wise  extent)  relative  to  flowfield  'noise',  an  axisymmetric  con¬ 
figuration  shown  in  Figure  1.  was  specified.  A  right  circular  cylinder,  whose  axis  was 
aligned  with  the  on-coming  flow  and  around  which  it  could  be  rotated  was  chosen,  in  part, 
due  to  its  relevance  to  various  U.S.  Air  Force  research  initiatives. 

Along  lhe  control  volume  surfaces  which  contained  the  cylinder,  sufficient  numerical 
boundary  conditions  were  specified  to  obtain  a  unique  solution.  At  a  large  radial  distance 
from  the  cylinder  uniform  freestream  flow  was  imposed.  At  the  downstream  end  of  the 
computational  domain,  outflow  boundary  conditions  were  used  to  simulate  a  cylinder  of 
at  least  semi-infinite  axial  extent.  With  respect  to  the  leading  'edge',  a  considerable 
number  of  test  cases  were  studied.  Initially,  various  specific  forebodies  were  employed. 

As  a  result,  the  flowfield  was  dominated  by  bow  shocks  and  expansion  waves,  bogging 
dcwn  the  numerical  analysis  on  interesting  phenomena  which  were  irrelevant  to  the 
researcn.  After  several  weeks  of  experimentation,  it  was  decided  to  abandon  the  approach 
of  employing  any  forebody  at  all.  It  was  assumed  that  the  cylinder  was,  in  fact,  infinitely 
long.  However,  in  the  region  -  *  <  x  <  0,  the  flow  was  allowed  to  slip  along  the 


surface  whereas  the  surface  temperature  equalled  the  freestream  static  temperature. 

In  this  manner  a  uniform  upstream  flow  entering  the  computational  domain  was  esta¬ 
blished.  In  the  region  0  <  x  <  10  Ax,  the  (axial  and  azimuthal)  velocity  slip  was 
linearly  decreased  to  zero  so  as  to  initiate  the  growth  of  the  boundary  layer  without  the 
very  strong  leading  edge  disturbances,  otherwise  computed  at  high  speeds,  that  previously 
dominated  the  downstream  flow. 

» 

With  this  set  cf  equations,  numerical  integration  scheme,  configuration  and  boundary 
conditions,  the  flowfield  was  computed  over  a  range  of  Mach  numbers  in  order  to  esta- 
^  blish  the  feasibility  of  utilizing  direct  numerical  simulation  for  computing  the  essential 

turbulent  shear  flow  structures/physics  unique  to  hypersonic  flow. 

» 

* 

b 
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Over  a  range  of  Mach  numbers  from  1.4  to  6.3  the  aforementioned  flowfield  was  com¬ 
puted  at  Reynolds  numbers  which  were  known  to  be  unstable  and  provide  significant 
viscous  interaction,  with  a  wail  temperature  equal  to  the  freestream  static  temperature. 
At  a  Mach  number  of  6.3  numerical  instabilities  could  not  be  controlled  with  a  time  step 
of  practical  size.  In  addition,  for  all  of  the  remaining  Mach  numbers,  the  low  numerical 
'noise'  resulted  in  unacceptably  long  computational  times  for  the  small  ambient  distur¬ 
bances  to  follow  linear  instability  growth  into  the  non-linear  three  dimensional  regime. 
During  pnase  II  various  ambient  noise  models  shall  be  added  to  the  simulation.  Therefore, 
all  of  the  numerical  simulations  presented  here  are  associated  with  an  artificial  excita¬ 
tion  described  earlier  in  the  form  of  a  temporally  harmonic  normal  wall  mass  transfer, 
with  periodic  suction/biowing  and  various  frequencies  relative  v'  the  time  scales  of  the 
mean  flow  and  viscous  layer  thickness. 

The  main  body  of  the  simulations  in  Figures  2-6  &  9-11  pertain  to  a  non-rotating  cylin¬ 
der.  Two  cases  including  cylinder  rotation  at  a  Mach  number  of  3.0  are  presented  in 
Figures  7  &  8  to  demonstrate  the  capability  for  including  three  dimensional  mean  flow 
effects.  Inasmuch  as  the  cylinder  rotation  results  did  not  qualitatively  affect  the  princi¬ 
pal  discovery  of  this  study,  they  shall  not  be  referred  to  further  and  are  left  for  the 
reader  to  review,  with  regard  to  the  non-rotating  cylnder,  the  results  are  presented  in 
groups  for  each  Mach  number  and  excitation  frequency  combination,  beginning  with  the 
highest  Mach  number,  within  each  set  of  simulations  for  a  given  Mach  number,  the  exci¬ 
tation  frequencies  include  one  or  both  of  the  cases  10"  1  and  10 1  times  the  frequency 
computed  from  the  local  freestream  velocity  and  boundary  layer  thickness.  These  are 
referred  to  as  low  and  high  frequency,  respectively.  The  cases  with  an  excitation 


frequency  cf  1QC  times  the  reference  were  essentially  identical  with  1CT!.  Within  each 
group  of  resuits  for  a  given  Mach  number  and  excitation  frequency,  the  fluctuating 
entropy,  azimuthal  component  of  vorticity  and  pressure  are  presented  along  with  charac¬ 
terizations  of  the  mean  flow.  Some  fifty  boundary  layer  thicknesses  downstream  of  the 
excitation,  the  effect  of  such  disturbances  on  the  flow  was  'observed*  and  is  presented 
in  the  various  plots. 

Although  it  is  true  that  the  bulk  of  the  research  time  was  consumed  with  structuring 
the  computations,  a  significant  observation  has  been  made.  This  is  most  readily  seen 
by  comparing  Figures  2a  &  2b  with  7a  &  7b.  As  the  Mach  number  increases,  the  entropy 
fluctuations  exhibit  a  second  distinct  boundary  layer  'edge'  well  beyond  that  established 
by  the  vorticity.  At  lower  Mach  numbers  the  entropy  fluctuations  radiate  much  further 
into  the  freestream  flow.  (As  a  part  of  Phase  II,  the  full  range  of  high  subsonic  and  tran¬ 
sonic  Mach  numbers  shall  also  be  evaluated  to  understand  the  entropy  fluctuation  field 
in  tne  regime  where  entropy  and  pressure  fluctuations  are  just  beginning  to  be  signifi¬ 
cant.)  As  a  result,  Phase  I  of  this  SB!R  program  has  not  only  demonstrated  the  feasibility 
of  utilizing  direct  numerical  simulations  to  expose  the  essential  physics  of  turbulent 
structure  in  hypersonic  flow,  but  has  a'so  provided  an  initial  direction  for  future  research 
into  the  possible  multiple  layer  structure  of  hypersonic  TB1_  fluctuations. 


Conclusions  and  Recommendations 


Numerical  simulations  of  boundary  layer  flow  from  the  hypersonic  down  through  tran- 

)  sonic  regimes  have  been  carried  out  utilizing  the  full  time-dependent,  viscous  and 

compressible  flow  Navier-Stokes  equations  for  a  perfect  gas.  The  nominally  axial  flow- 
field  over  a  circular  cylinder,  with  and  without  rotation  about  its  axis,  was  computed 

i  utilizing  existing  numerical  techniques  coupled  with  a  unique  leading  edge  flow  develop¬ 

ment  formulation.  In  order  to  generate  a  numerically  'quiet'  solution,  the  (axial  and 
azimuthal)  zero  slip  boundary  conditions  were  introduced  gradually  over  the  first  ten 

>  grid  points  beginning  at  the  leading  edge  of  the  cylinder.  A  linear  distribution  ranging 
from  full  slip  at  the  first  node  to  zero  slip  at  the  tenth  node  was  imposed.  The  wall 
temperature  was  set  equal  to  the  freestream  static  temperature.  In  future  simulations 

>  with  various  hot  and  cold  walls  the  temperature  difference  shall  also  be  phased  in  over 
the  first  ten  nodes.  This  procedure  allows  for  the  computation  of  a  generic  boundary 
layer  flow  which  is  essentially  noise  free,  without  the  complications  associated  with  the 

)  specification  of  a  particular  forebody. 

Due  to  the  low  noise  level  in  the  flow,  growth  of  natural  disturbances  into  instability 

>  waves  with  subsequent  non-linear  three-dimensional  breakdown  required  intolerably  long 
computational  times.  In  order  to  meet  the  Phase  I  SBIR  cost  and  time  constraints,  the 
boundary  layers  were,  therefore,  artificially  excited  with  a  temporally  harmonic  distur- 

>  bance  at  the  wall  several  boundary  layer  diameters  upstream  of  the  observation  station. 
Both  periodic  suction  and  blowing  as  well  as  heating  and  cooling  were  evaluated,  with 
the  former  being  selected  for  this  study.  The  peak  magnitude  of  the  mass  flux  was  held 

>  constant  at  1%  of  the  freestream  whereas  the  time  constant  was  varied  over  a  range  of 
values  relative  to  the  viscous  flow  length  and  velocity  scales. 

> 


By  utilizing  sucn  artificial  excitation,  the  development  of  disturbances  into  the  non¬ 
linear  regime  and  the  beginnings  of  essentia!  three-dimensionality  could  be  computed 
over  a  wide  range  of  Mach  numbers.  For  the  sake  of  economy  of  time  and  costs  the 
computations  were  terminated  when  any  essential  characteristic  of  the  disturbance  field 
was  observed  at  tnat  Mach  number  so  that  a  wide  range  of  Mach  numbers,  frequencies 
and  mean  flow  conditions  could  be  reported  upon. 

As  a  result  of  these  computations,  a  unique  effect  of  Mach  number  on  flowfield  distur¬ 
bances  has  been  observed,  when  viewed  in  terms  of  the  vorticity,  entropy  and  pressure 
fluctuations,  hypersonic  flows  appear  to  exhibit  a  distinct  multiple  layer  structure  in 
which  the  vorticity  fluctuations  are  confined  to  the  region  established  by  the  mean 
vorticity  whereas  the  entropy  and  pressure  fluctuations  extend  further  out  into  the  uni¬ 
form  flow  yet  still  terminate  relatively  near  the  solid  boundary.  At  lower  supersonic 
Mach  numbers  the  entropy  and  pressure  fluctuations  propagate  increasingly  far  away 
from  the  body  and  are  reminiscent  of  shadowgraph  pictures  of  waves  emanating  from 
the  eoge  of  a  turbulent  boundary  layer. 

In  accordance  with  the  objectives  of  a  Phase  I  5BIR,  the  feasibility  of  utilizing  direct 
numerical  simulations  to  expose  the  essential  physics  of  turbulent  boundary  layer  struc¬ 
tures  at  hypersonic  speeds  has  been  demonstrated.  It  is  apparent  that  simulations  carried 
out  in  this  manner  could  be  invaluable  in  both  the  complete  description  of  the  unique 
vcrticity-entropy-pressure  field  dynamics  as  well  as  providing  guidance  in  the  formula¬ 
tion  of  more  compact  engineering  models. 

Based  on  this  established  feasibility,  an  exhaustive  computational  study  is  indicated  in 
which  longer  computational  times  could  be  tolerated  and  natural  disturbances  followed 
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into  non-linearity  and  eventually  chaos.  It  would  be  particularly  advantageous,  if  the 
Phase  II  SBIR  effort  is  able  to  be  supported,  if  the  work  would  be  associated  with  an 
effort  in  which  experimental  data  on  hypersonic  turbulent  flow  is  being  obtained.  The 
Hokenson  Company  is  under  contract  to  AFATL/MNG,  Eglin  AFB,  FL  to  carried  out 
experimental  studies  on  spinning  hypersonic  projectiles  with  base  bleed  flow  and  optimal 
forebody  shaping,  it  is  our  recommendation  that  the  Phase  II  work  on  computing  turbu¬ 
lent  hypersonic  flow  boundary  layer  and  wake  structures  be  linked  to  the  geometrical 
configuration  being  studied  therein,  in  this  manner  two  research  programs  which  are 
exclusively  computational  and  experimental,  respectively,  would  benefit  from  the  associ¬ 
ation  with  one  another  and  the  achievements  of  two  Air  Force  research  interests  would 
be  significantly  enhanced. 
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^  Supplemental  Information 

The  professional  personnel  involved  with  the  work  were  Dr.  Gustave  J.  Hokenson  and 
^  Mr.  Richard  D.  Crouse,  who  holds  an  MS  in  Aerospace  Engineering  from  Stanford 

University. 

9  Contingent  upon  obtaining  Phase  I!  SBIR  support,  The  Hokenson  Company  intends  to 

present  several  papers  related  to  this  work  at  the  First  National  Fluid  Dynamics 
Congress  with  subsequent  journal  publication. 

b 


Results  of  this  effort  have  also  been  invaluable  in  supporting  an  SBIR  program  being 
carried  out  by  The  Hokenson  Company  for  AFAtl/mng,  Egiin  AFB.  In  that  work,  drag 
reduction  on  20- mm  hypersonic  projectiles  is  being  studied  experimentally  through  appli¬ 
cation  of  base  bleed  and  optimal  foreoody  shaping.  Advances  in  the  understanding  of 
turbulent  transport  processes  derived  from  the  AFOSR  Phase  1  SBIR  work  reported  upon 
here  significantly  impacts  the  design  of  an  effective  experimental  program.  ]n  addition, 
the  experimental  data  from  such  a  study  would  significantly  enhance  the  ability  to  vali¬ 
date  the  numerical  simulations  scheduled  to  be  carried  out  during  the  AFOSR  Phase  II 
SBIR  effort. 
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Case  I.:  M=5.3,  Zero  Rotation 


Fiqure  No. 


2a.  Entropy  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation. 
2d.  Vorticity  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Escitation. 
2c.  Pressure  PJuctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation. 
2d.  Mean  Flow  Mach  Number  Distribution . 


Entropy  Fluctuations  vs.  Radius  and  Time,  High  Frequency  Excitation, 
vorticity  rluctuations  vs.  Racius  and  Time,  High  Frequency  Excitation. 
Pressure  Fluctuations  vs.  Radius  and  Time,  High  Frequency  Excitation. 
Mean  Plow  Mach  Number  Distribution . 
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Entropy  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation.  .  .  39 

Vorticity  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation.  .  .  40 
Pressure  fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation.  .  .41 
Mean  Flow  Macn  Number  Distribution . 4  2 
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Case  Ilia.;  M=3.0,  Zero  Rotation 


p'iaure  No. 


3a.  Entrooy  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation. 
5b.  Vorticity  fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation 

5c.  Pressure  fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation. 

5d.  Mean  Flew  Mach  Number  Distribution . 

6a.  Entropy  Fluctuations  vs.  Radius  and  Time,  High  Frequency  Excitation. 
6b.  Vorticity  Fluctuations  vs.  Radius  and  Time,  High  Frequency  Excitation 

6c.  Pressure  Fluctuations  vs.  Radius  and  Time,  High  Frequency  Excitation 

6d.  Mean  Flow  Mach  Number  Distribution . 
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Case  I!Ib.:  M=3.0,  Non-Zero  Rotation 


Figure  No. 


Entropy  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation,  Surface 

Rotational  Speed  =  .25  Freestream . 5 

Vorticity  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation,  Surface 

Rotational  Speed  =  .25  Freestream . 5 

Pressure  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation,  Surface 

Rotational  Speed  =  .25  Freestream . 5 

Mean  Flow  Swirl  Distribution . 5 

Entropy  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation,  Surface 

Rotational  Speed  =  .5  Freestream . 5 

Vorticity  fluctuations  vs.  Radius  and  Time,  Low  Frequency  Ex  citation,  Surface 

Rotational  Speed  =  .5  Freestream . 5 

Pressure  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation,  Surface 

Rotational  Speed  =  .5  Freestream . 5 

Mean  Flow  Swirl  Distribution . 6 
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Case  IV.:  M=2.Q,  Zero  Rotation 


Figure  No. 

9a.  Entropy  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation. 
90.  Vcrticity  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation 
9c.  Pressure  Fluctuaitons  vs.  Radius  and  Time,  Low  Frequency  Excitation 
9d.  Mean  Flow  Mach  Number  Distribution . 


.“in.' 


nson  Company 


enson  Company 


o  o 


Case  V.:  M-1.4,  Zero  Rotation 


‘iaure  No. 


ua.  Entropy  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation, 
do.  vorticity  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation 

c.  Pressure  Fluctuations  vs.  Radius  and  Time,  Low  Frequency  Excitation. 

d.  Mean  Flow  Macn  NjnnOer  Distribution . . 

ia.  Entropy  Fluctuations  vs.  Radius  and  Time,  High  Frequency  Excitation. 
id.  vorticity  Fluctuations  vs.  Radius  and  Time,  High  Frequency  Excitation 
lc.  Pressure  Fluctuations  vs.  Radius  and  Time,  High  Frequency  Excitation 
id.  Mean  flow  Mach  Numoer  Distribution . . 
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